Cosmological magnetic fields from nonlinear effects 
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In the standard cosmological model, magnetic fields and vorticity are generated during the ra- 
diation era via second-order density perturbations. In order to clarify the complicated physics of 
this second-order magnetogenesis, we use a covariant approach and present the electromagneto- 
dynamical equations in the nonlinear regime. We use the tight-coupling approximation to analyze 
Thomson and Coulomb scattering. At the zero-order limit of exact tight-coupling, we show that 
the vorticity is zero and no magnetogenesis takes place at any nonlinear order. We show that 
magnetogenesis also fails at all orders if either protons or electrons have the same velocity as the 
radiation, and momentum transfer is neglected. Then we prove a key no-go result: at first-order in 
the tight-coupling approximation, magnetic fields and vorticity still cannot be generated even via 
nonlinear effects. The tight-coupling approximation must be broken at first order, for the generation 
of vorticity and magnetic fields, and we derive a closed set of nonlinear evolution equations that 
governs this generation. We estimate that the amplitude of the magnetic field at recombination on 
the horizon scale is ~ 1Q -27 G. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The origin of cosmological magnetic fields is an impor- 
tant problem in cosmology [1]. Many mechanisms for 
primordial magnetogenesis (i.e., creation before struc- 
ture formation) have been proposed. In order to gen- 
erate fields on large scales, inflationary mechanisms are 
the best candidates, but they require uncertain modifica- 
tions to standard physics in order to break the conformal 
invariance of Maxwell fields [2] . 

The generation of cosmological magnetic fields via 
plasma interactions during the radiation era, originally 
suggested by Harrison [3], is based on conventional 
physics and does not require any new postulates. The 
essential ingredients in this mechanism are nonzero vor- 
ticity and Thomson scattering between photons and 
charged particles. Because momentum transfer is more 
effective between photons and electrons than between 
photons and protons due to the mass difference, Thom- 
son scattering induces differences in the velocity and the 
distribution of protons and electrons. These differences 
induce local electric currents and net charge density, and 
the electric field in turn generates a magnetic field. This 
process initiates after electron-positron annihilation and 
ends when there are insufficient free electrons, i.e., it op- 
erates over the temperature range T roc < T < m e . 

In the standard cosmological model, there is no vortic- 
ity at first order, and the generalized Harrison mechanism 
is much more complicated. We can follow the basic argu- 
ment in a simple Newtonian formalism. The evolution of 
the magnetic field is described by the induction equation, 

B = -VxE. (1) 
Analysis of the momentum transfer in scattering leads to 



the generalized Ohm's law, 

E = T)J + §, 



(2) 



where r\ is the plasma resistivity and S is the contribu- 
tion from Thomson scattering. Since r\ <C H -1 , we can 
neglect this term to obtain 



B = —V x S . 



(3) 



As we will see below, S ~ nv, where n is the number den- 
sity of charged particles and v is the velocity difference 
between radiation and charged particles. Equation (3) 
shows that there are two sources for magnetogenesis - 
vorticity, curl v, and the vector product of density gradi- 
ent and velocity, Vn x v. 

In the standard cosmology, where perturbations are 
generated from inflation, there are no vector modes at 
first-order, and therefore the vorticity vanishes at first 
order. The density-velocity term is a product of first- 
order scalar perturbations and therefore also vanishes at 
first order. Thus in the standard model a perturbative 
analysis of magnetogenesis during the radiation era must 
start at second order. (Exceptions can arise if there are 
sources of first-order vector perturbations, such as cos- 
mic strings [4], or fine-tuned anisotropics in collisionless 
neutrinos [5].) 

Matarrese et al. [6] analysed how vorticity and mag- 
netic fields can be generated from second order cosmo- 
logical perturbations. Subsequent work has also used 
second-order perturbations [7-11], but has neglected the 
vorticity and metric vector perturbations, and focused 
on the density- velocity terms, i.e., the product of first- 
order scalar terms. (For other work on magnetogenesis 
during the radiation era, see Refs. [12].) The different 
approaches lead to estimates in the range [6, 7, 10, 11]: 



B ~ 10 



-24 



10~ 27 G 



(4) 
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at recombination on 100 Mpc scales. This is a very weak 
field, but it provides a seed which is amplified via the 
dynamo mechanism. It is possible that the dynamo am- 
plification can reach the current observed value of about 
10 _6 G on galaxy scales [13]. 

The simplistic Newtonian description given above al- 
lows us to identify the key physical effects, but the real 
situation is much more complicated. The second-order 
perturbative treatments arc a necessary foundation for 
computing the power spectrum of the magnetic field. 
However, it is also useful to adopt a covariant approach 
that directly generalises the Newtonian treatment to cos- 
mology [14, 15]. This allows us to develop a direct 
physical understanding of the magnctogenesis process, 
and also to deal with the problem in the fully nonlin- 
ear regime. 

The greatest complexity arises from the dynamics of 
momentum transfer via scattering. We use the tight cou- 
pling approximation [16], which is based on the fact that 
the scattering time r is much less than the cosmic expan- 
sion time if -1 , 

Ht < 1 , (5) 

so that photons and charged particles are closely bound. 
In the limit, i.e., at the zero-order of exact tight-coupling, 
we have r = and v = 0, so that all particles share the 
same velocity and behave as a single fluid, and no mag- 
netogenesis takes place. Beyond the zero-order of tight 
coupling, there is a nonzero t and velocity difference v, 
and the tight coupling approximation is an expansion in 
Ht. Note that the tight coupling approximation is in- 
dependent of the cosmological perturbative approxima- 
tion. Zero-order in tight coupling is not to be confused 
with zero-order in cosmological perturbations; the cos- 
mological variables can be at any nonlinear order, but no 
magnetogenesis is possible. 

First, we give a basic result: in the limit of exact 
tight coupling, and neglecting anisotropic stresses, vor- 
ticity vanishes and magnetic fields cannot be generated 
at any nonlinear perturbative order. If the exact tight 
coupling limit is partially relaxed by assuming that ei- 
ther the protons or the electrons have the same veloc- 
ity as the radiation, but neglecting momentum transfer, 
then vorticity and magnetic fields still cannot be gener- 
ated at any nonlinear order. Then we derive the evolu- 
tion equation for magnetic fields and vorticity beyond the 
zero-order of the tight coupling approximation. We show 
that there is no magnetogenesis at the first order in the 
tight- coupling approximation, and magnetogenesis takes 
place at second order. While it has been suggested that 
magnetogenesis is possible via the breaking of the tight 
coupling limit [17], this is the first study giving explicitly 
the condition for such a mechanism to work. 

The rest of the paper is organised as follows. In the 
next section we derive the nonlinear equations for mag- 
netic fields and vorticity in a covariant formalism. In 
section III, we show that vorticity cannot be generated 
even through nonlinear effects in the tight coupling limit. 



In section IV we consider the first and second order of 
the tight coupling approximation. Finally we summarise 
our work in section V. 



II. COSMOLOGICAL 
ELECTROMAGNETO-DYNAMICS 

The Faraday tensor can be split into electric and mag- 
netic fields as measured by a congruence of fundamental 
observers u a (with u a u a = —1): 

F a b — 2u[ a E b ] + e a b c B c , (6) 

where E a u a — B a u a = 0. The spatial alternating tensor 
is s a bc = i] a bcdU d , where r\ abc d is the spacetime alternating 
tensor, using the convention 770123 = —\f-~9- The tensor 
indices represent an arbitrary coordinate or tetrad frame; 
at any event one can choose local inertial coordinates 
such that u a = (1, 0), E° = = B°. 

The induced metric in the observer's comoving rest 
space is 

h a b = 9ab + u a u b , (7) 

and it defines a covariant spatial derivative D a . Gener- 
alizing the Newtonian case, we define kinematical quan- 
tities of the u a congruence via its covariant derivative: 

V h u a = ^8h ab + a ab + e abc uj c - u a u b . (8) 

Here 9 = V a u a is the volume expansion rate, a ab — 
h a c h b d \7 ( c u d ) is the shear, uo a = — ^e abc \7 b u c is the vor- 
ticity, and ii a — u b V b u a is the acceleration. 

Since we are working mainly in the radiation era, it is 
useful to choose u a as the radiation four-velocity in the 
energy frame, i.e., with no energy flux: 

u a = u a 7 , q" := -T^ b u b - p^u* = . (9) 

The four-velocities of charged particles, I — p,e, are 

«? = 7/ k + vn , Ua Vj a - , 7/ = (1 - v Ia vn~ 1/2 , 

(10) 

where we also choose the energy frames, so that qf := 
-Tf b u b -p I u a I = Q. 

The Maxwell equations V [a F fcc] = and V b F ab = j a 
can be split in a 1+3-covariant way (relative to u a ) as [14] 

D a B a = 2E a w\ (11) 
B a E a = -2B a u a + ti, (12) 

B a = -^9B a + (a ab + e a bcio c )B b 

- curl E a -e abc u b E c , (13) 

E a - ~6E a + {o- ab + e abc LJ c )E b 

+ curl B a + s abc u b B c - J a . (14) 
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Here n = —j a u a is the charge density and J a = h ab j b is 
the current. For the radiation era plasma with T < m e , 

f = Jp+Je, j? = e I n I u a I , (15) 
fj, = e(7 p n p - j e n e ) , (16) 
J a = e( lp n p V p a - le n e V e a ), (17) 

where e\ = ±e, and nj are the number densities. We can 
write 



m = n(l + Aj), 



(18) 



where n is the density of charged particles in the tight 
coupling limit. 

The four-current satisfies local charge conservation, 
V a j a = 0, which implies 



(1 + 0,1 = -B a J a -u a J a . 



(19) 



In order to close the Maxwell equations, we need to spec- 
ify j a , and this is done via the equations of motion for 
photons and charged particles. 

In the Maxwell equations we use the covariant curl, 



curlS a :=e abc D b S c , 



(20) 



and the ovcrdot is the covariant time derivative along the 
radiation 4- velocity u a , projected into the rest-space: 



S a := h b u c V c S b . 



(21) 



(This is more convenient than the usual definition, which 
is not projected.) At any spacetime event one can 
choose inertial coordinates so that D a / = (0,Wf)S a = 
(0,d t S),cur\S a = (0, V x S). Two important identities 
are [15] 

D a / = (D„/)" -fu a 

+ (\0h ab + a ab - e abc u c \ D 6 /, (22) 



curlD a / 



-2/w„ • 



(23) 



The vorticity propagation equation is independent of 
the field equations, and is given covariantly by [15] 



2 . 1 

UJ a = -g^a + CTafcW - -CUrl?i a . 



(24) 



The conservation law for the electromagnetic energy- 
momentum tensor Tf = F a c F bc - \F cd F cd g ab follows 
from the Maxwell equations: 

V b T£ b = -F ab Jb . (25) 

Photons and charged particles obey the balance equa- 
tions 



V b T 7 aft = K« , VbTf = Kf + e ini F a b u b , (26) 

where the K a four-vectors are the rates of energy- 
momentum density transfer to the species. By 



Eqs. (15) and (25), the conservation of the total energy- 
momentum, \7 b (Tf + Y,i T f + t< f) = °> implies that 
K« + K$ + K% = 0. 

The photon energy and momentum balance equations 
in the general nonlinear case are [15] 



Pi + = 
- Pl u a + ^D a Pl = 



ab 

7 



II 



D b 7T° b — u b TrH. b + M® , (28) 



7 ' 



ab 



(27) 



where 7r" fc is the anisotropic stress, with n^ b u a = = 



rf b h ab . 



Here XJ 1 = -u a K* and 



h a b K b are 



the rates of energy and momentum density transfer to 
photons from Thomson scattering. From now on, we 



take 



T ab 



0; the role of photon anisotropic stress in 



magnetogenesis has been investigated by Takahashi et 
al. [8, 9, 11]. For electrons and protons, it is reasonable 
to neglect pressure and anisotropic stresses. Then the 
energy conservation equations are 



U/V 6 p/ + e lPl = Uj , 



(29) 



where pi = mjnj and Ui = —uj a Kf is the rate of energy 
density transfer due to Thomson and Coulomb scatter- 
ing. The momentum balance equations are 



p/U^V&u? = Mf + e ini F a b u b j , 



(30) 



where Mf = hq h K b is the rate of momentum density 
transfer due to Thomson and Coulomb scattering. 

As shown in Maartens et al. [15], the Thomson energy 
transfer rates, U 7 and Ui, start from 0(Vf) and 0{Vj), 
respectively, while the Thomson momentum transfer 
starts from linear order in Vj. As we will explain later, all 
0{Vf) terms, except those from Thomson scattering, can 
be neglected in order to derive our evolution equation for 
magnetic fields. Writing down the necessary terms ex- 
plicitly, we have 

U 7 = ]TC 7/ V} 2 + 0(C 7/ V} 3 ), (31) 

Ui = 0(C 7l Vj 3 ), (32) 
M i = -£<V(u 6 -«i)fc°6 

= ^C 7/ V 7 a + 0(C 77 y/), (33) 

Mf = -C lI {u b I -u b )h a Ib -C I j{u b I -u b J )h a Ib 
= -C^V? + V?u a ) -Cu{Vf -Vf) 

+0{C lI Vf) , (34) 

where C 7 j, Cjj are the Thomson and Coulomb collision 
coefficients. The energy and momentum balance equa- 
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tions reduce to 

Pi + \pi® = E C -t' V i + °( C H V i) > 
6 i 

h T + nj(e + D a V T a + u a V T a ) + VfDam 

= o(c 7/ y/) + o(v?) , 



- Pl u 



k a + Vp 7 =J2 C H V i a + 0(C 7 jV?) , 



(35) 

(36) 
(37) 



m ini (ii a + u b V b V? + V b V b u a ) 

= e ini (E a + F\Vj b ) - C 7l Vj a - Cu(Vj a - VJ 
+ 0(C 7l V?) . 



(38) 



The vorticities of charged particles are related to the ra- 
diation vorticity uj a at O(Vj) by [15] 

uj = uJ a - icurl Vf + ^e a bc u b Vf + u b Vfu a . (39) 

Taking the tight coupling limit in Eq. (36) we recover 
number conservation, 



n + 6n = . 



(40) 



and with Eq. (35), this leads to w a V a [ln(n//» 
We define the entropy 



Sa ■= D a 



Using the identity (22), we arrive at 



s a + ^0s a + {o ba + e bac uj c )s b = . 



3 / 4 ^l - 




= 0. 



(41) 



(42) 



In what follows we assume the adiabatic condition s a = 0, 
i.e., 



D a n _ 3 D a p 7 
n 4 p 7 ' 

which is consistent with Eq. (42). 

III. TIGHT COUPLING LIMIT 



(43) 



The induction equation (13) reduces to 



B a — --8B a + a ah B h . 



(45) 



so there is no source term and no magnetogenesis. Equa- 
tions (44) and (45) hold in the fully nonlinear regime. 

Next we consider what happens if exact tight coupling 
is weakened by neglecting scattering terms, and neglect- 
ing the velocity difference between protons and photons, 
i.e., U p a = 0, but allowing U e a ^ 0. This is effectively the 
assumption made in Ref. [6], and here we reconsider the 
problem in the covariant formalism. The proton equation 
of motion (38) becomes m p u a = eE a , and taking the curl 
gives 



curlu a = curl E a 



(46) 



The curl of the photon momentum equation (37), using 
the photon energy equation (35) and the identity (23), 
gives 



curlu a 



(47) 



Using these equations, the induction equation (13) and 
the vorticity propagation equation (24) become 

2 1m 

K = -T0B a + {a ab + e abc u c )B b + — ^6u a , (48) 

3 3e 



(49) 



Thus vorticity and the magnetic field are conserved and 
no magnetogenesis is possible. This gives our first no-go 
result: 

If (i) the scattering terms are neglected and the proton- 
photon velocity difference is neglected, (ii) anisotropic 
stress is neglected, (Hi) the initial magnetic field and vor- 
ticity are zero, and (iv) the energy-momentum conserva- 
tion equations hold, then vorticity and magnetic fields 
cannot be generated, at any perturbative order. 

The same result holds if we assume the electron veloc- 
ity equals the radiation velocity but instead Vp ^ 0. It 
is not clear how the results in Ref. [6] relate to our no-go 
result. 



In the exact tight coupling limit, the velocities of pho- 
tons, protons and electrons are equal, Vf = 0, and the 
momentum transfer terms vanish. Intuitively we expect 
that vorticity and magnetic fields cannot be generated 
from zero. In fact, we can prove a stronger result, that 
vorticity is zero at all nonlinear orders. Equation (38) 
reduces to mjii a = eiE a , which implies E a = = u a . 
Then the photon momentum balance equation (37) re- 
duces to D a /9 7 = 0, and since p 7 ^ by Eq. (35), the 
identity (23) implies 



IV. TIGHT COUPLING APPROXIMATION 

Nonlinear magnetogenesis is ruled out in the tight cou- 
pling limit of zero collision time, r = 0. Beyond the 
zero-order of tight coupling, there is a nonzero r and a 
nonzero velocity difference v, which is governed by the 
momentum balance equations (37) and (38). Schemati- 
cally, these are of the form 



U a = 0. 



(44) 



v = - + A, 

T 



(50) 
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where A represents terms other than scattering terms. 
Since v <~ Hv, we have v <C vjr. We expand in terms of 
the tight coupling parameter tH: 



v = v {1) + v {2) H , A = A {0) + A {1) 



(51) 



and we use TCA(n) to denote n-th order in the tight 
coupling approximation. Then 



TCA(l): 0=^ + A (0) , 
TCA(2): 6 (1) = M+^ (1) . 



(52) 
(53) 



The TCA is complicated by the presence of Coulomb 
scattering, so that strictly we need to perform TCA ex- 
pansions in both Thomson and Coulomb small parame- 
ters. However, as we will argue, it is reasonable to neglect 
the Coulomb collision time, i.e., to assume tight coupling 
of protons and electrons. 

The collision coefficients in Eqs. (37) and (38) are 



C 7 e — gOTP 7 n e , 

r - d 2 ^-C 

7P — ji ' e ' 

C pe = e 2 n e n p r) , 
where <7t is the Thomson cross-section and 



(5:= 



rn 



p 



Here rj is the resistivity of the cosmic plasma, 



V : 



4ire 2 /m e \ 3 / 2 



V T 



In A - 10~ 13 | — 



in, L ) " VeV 



-3/2 



(54) 
(55) 
(56) 

(57) 
(58) 



where In A is the Coulomb logarithm, A ~ T 3 / 2 e _3 ri _1 / 2 
and n ~ 10~ 10 T 3 . The C's define key timescales, together 
with the Hubble timcscale: 



T<VT> • 



' ep •- 



m e n e 
r 

m e n e 

,12 



(59) 



T 



- -10 ia ( — ) s, (60) 



^'eV 



T 



-3/2 



Thus 



Hr ie ~ 10- 7 ( - 



1/2 



(61) 
(62) 

(63) 
(64) 
(65) 



As one can see, Thomson scattering between photons and 
electrons and Coulomb scattering between protons and 
electrons are very effective on cosmological timescales, 
so that they are tightly coupled before recombination. 
Although Thomson scattering between photons and pro- 
tons is less effective at low temperatures, protons also 
closely follow photons through their Coulomb coupling 
to electrons. 

From Eqs. (64) and (65), we see that, 



^£ ~ 10 9 f — 

ev 



' ep 



-5/2 



(66) 



Therefore, at lower temperatures, T < 1 keV, Coulomb 
scattering is more effective than Thomson scattering, so 
that protons and electrons are more tightly coupled than 
photons and charged particles are. This suggests that we 
can safely neglect Coulomb scattering, i.e., take r ep = 







ya 
"p 



V, a = A v — A e . We define the centre of mass 



velocity V a and number density deviation A: 



(m p n p + m e n e )V a = m p n p V p a + m e n e V^ , (67) 
(m p + m e )A = m p A p + m e A e . (68) 

Then the peculiar velocities decompose as 



ya = ya + v a _ ya 



(69) 



where u" are the deviations of proton and electron ve- 
locity from their centre of mass velocity, with m p n p v p + 
m e n e v^ = 0. The number density deviations decompose 
as 



A/ = A + (5 7 w A, 



(70) 



where 5i are the deviations of proton and electron num- 
ber density from their centre of mass density, with 
m p 5 p + m e 5 e = 0. The approximations in Eqs. (69) and 
(70) are based on neglecting terms of order r ep : 

iJr 7e ~ |A| ~ \V a \ » H Tep ~ 1(5/1 - K| . (71) 

In fact, protons and electrons are coupled not only by 
Coulomb scattering but also by the electric field, so that 
\Si\, \v"\ are further suppressed by a factor [18], 



Hr) - 10 



-27 



T 
cV 



1/2 



(72) 



Thus we can safely apply the approximation at any tem- 
peratures we consider here, m e > T > T rcc , and we can 
solve equations perturbatively with respect to the small 
parameter f/T 7e . 

In Eqs. (36)-(38), we can set V? = V a . We can see 
from Eqs. (52) and (53) that, upto TCA(2) we consider 
here, 0(C 7 /Vj 3 ) and 0(V 2 ) terms can be neglected and 
the momentum equations become, 

^>" + l D > 7 = (/?3m ^ + mene V », (73) 



m 7 (ii a + u b \7 b V a + V h \7 b u a ) 



= ej (E a + F a b V h 



mi 
ti 



V a . 



(74) 
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where we have defined 

T := T 7e , 77 := T-yl = (/3~ 3 T, t) . 



(75) 



On the other hand, as we will see in section IV B, we 
need only TCA(l) for the charged particle conservation 
equation (36) to derive the evolution equation: 

A + B a V a +u a V a + V a — = 0. (76) 



A. First order - TCA(l) 

At first-order TCA, we follow Eq. (52) and keep only 
the first-order V a and the zero-order of other terms. 

By summing the photon and charged particle equa- 
tions (73) and (74), we obtain 



4p 7 (l + R) 



p 1 . 



where 



R := - — , p b := (m p + m e )n . 
4p~ f 

The charged particle equations (74) imply 




(i + 13)— £ a + (1 ^ v: 
m e t 



(i) 



(77) 



(78) 



(79) 



The photon and charged particle equations (73) and (74) 
also lead to 







l + P 2 \ 0(1 + R) 
1 + /3 



(80) 



It follows that 



_ m e (1 - (3 3 ) 1 

m p (1 - /3 3 ) 1 



e (1 + /J 2 ) 4p 7 (l + i?) 



D°p 7 . (81) 



Using Eqs. (77) and (81) in the induction equation (13), 
we find the evolution equation for the magnetic field at 
TCA(l): 



B a = --6B a + (o ab + e abc u c )B b 



m p (1 - /3 3 ) 



Pi 



e (1 + /3 2 ) 2 Pl {l + R) 



(82) 



where we used the adiabatic condition (43). The mag- 
netic field is sourced by the vorticity of photons. How- 
ever, from Eqs. (43) and (77), we see that curlw a = 
[p 7 /2p 7 (l + R)]uj a , and the vorticity propagation equa- 
tion (24) becomes 



2R p 7 



4(1 + R) Pl 



(83) 



This shows that there is no source term for the vorticity, 
and we have our second no-go result: 
No generation of vorticity or magnetic fields is possible 
at first order in the tight coupling approximation. 

Note that at TCA(l), terms of the form Vn x v do not 
arise. These terms come in at TCA(2). Magnetogene- 
sis via the V x v vorticity term is not possible without 
breaking the tight coupling approximation at first order; 
it is not clear how the results of Ref. [6] conform to this 
no-go result. 



B. Second order - TCA(2) 

Now we proceed to the second-order TCA version of 
the nonlinear evolution equations for the magnetic field 
and vorticity. 

The charged particle equations (74) imply 

E a + F a b yb = _H^ (1 ~ g lya (g4) 

e (1 + p) t 

Substituting into the induction equation (13), we have 



B a = --6B a + (a ab + s abc Lu c )B b 



+ 



m e (1-/3 3 ) 1 
e (1+/3) r 



cm l V "-Z e °bc Vb *±±L 
4 p 7 



,(85) 



where F a b V b in Eq. (84) was dropped because it is always 
negligible compared to B a . 

We need to solve for V a at TCA(2). The photon and 
charged particle equations (73) and (74) give 



u b V b V a + V b V b u a -\^^ 
4 Pl 

P(l + P 2 ) (1 + R + RA) ya 



(1 + /3) 

This can be split into TCA(l) and TCA(2) equations: 
1D> 7 _ /3(1 + /3 2 )(1 + J R) 



(86) 



4 p 7 



(i) ' 



(87) 



u h V h Vfa + VL\7 b u a 



'(i) -r- "(i) 
/3(1 + /3 2 ) 1 
(1 + /3) r 



(l + R)V ( %+RA {1) Vfa ,(88) 



where the first equation is equivalent to Eq. (80). We 
can solve these equations for V a : 



v a 



1 



(1) 



(1 + P) D> 7 

-T- 



4(1 + R) p(l + P 2 ) Pl 



(89) 



v In 



(2) 



R (1+/3) D> 7 
rA (1) - 



4(l + i?) 2 /3(l+/3 2 ) 



Pi 



1 



(1+P) 



(l + R) (3(1 + P 2 ) 



u^bV^+V^VbU 11 (90) 
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The last term on the right of Eq. (90) is determined in 
terms of D a p 7 from Eq. (89), and the vorticity occurs 
explicitly since VLV b u a oc B b Pl [\6h a b + a a h + e a bc u c ]. 
Now we can compute the crucial term in Eq. (85): 



curlU a 



£ abc Vb ^cPj 



Pi 

(1 + /3) 



2(1 + R) 13(1 + /3 2 ) p. ( 

R (1+/3) _abc 



D c/ 9 7 



4(1 + R) 2 13(1 + 



re at>c D b A (1) -^.(91) 
Pi 



Finally the evolution equation (85) for the magnetic 
field can be written, up to TCA(2), as 



B a = — 6B a + a a b B b 
3 



in 



v 



R 



(1-/3 3 ) 



e 4(1 + R) 2 (l + (3 2 ) 
m p 1 (1 - /3 3 ) p 7 



lbc D 6 A 



(i)- 



D c p 7 



Pi 



e 2(1 + R) (l + (3 2 ) Pl 



(92) 



The vorticity evolution equation (24) can be rewritten, 
using Eqs. (73), (89) and (90), as 



w = 



(1 + 2R) Pl 



A(l + R) Pl 
R 

8(1 + R) 2 



c D fc A (1) 



D c p 7 
Pi 



(93) 



The evolution of the baryonic number density deviation 
is governed by Eq. (76), which becomes, up to TCA(l), 



A <-'-> =-^(i)^- D ^(V< y (i)' ( 94 ) 



where Vfa is given by Eq. (89). 

Equations (92), (93) and (94), for given D a p 7 /p 7 , form 
a complete set of equations which describe the evolution 
of A(!), B a and uo a . We can see that both magnetic field 
and vorticity are generated at the second order in the tight 
coupling approximation. 



V. SUMMARY 

In this paper we have derived the evolution equations 
for cosmological magnetic fields and vorticity using the 
1+3-covariant formalism. The covariant approach al- 
lows us to construct a set of equations describing the 
fully nonlinear evolution of cosmic inhomogeneities. We 
have performed a tight coupling expansion for Thomson 
and Coulomb interactions to make the key physical pro- 
cesses transparent. It should be enphasized that we have 
not expanded inhomogeneous quantities with respect to 
cosmological perturbations, and therefore our results are 
valid at any order in cosmological perturbations. Thus, 



the present analysis is complementary to previous studies 
based on cosmological perturbation theory. 

Our first no-go result is that magnetic fields and vor- 
ticity cannot be generated in the tight coupling limit or 
its weak extension without anisotropic stresses. Then 
we have considered leading and next-to-leading order ef- 
fects in the tight coupling approximation. We have found 
that magnetic fields and vorticity are not generated at 
first order in the tight coupling approximation [TCA(l)]. 
The second order tight coupling approximation [TCA(2)] 
is necessary for generating both of them, and we have 
derived a closed set of nonlinear evolution equations at 
TCA(2). It is worth noting that we have not invoked the 
Einstein equations, so that our result does not rely on 
any specific theory of gravity. 

The magnitude of the generated magnetic field can be 
roughly estimated as follows. From Eqs. (89) and (94), 
we have A^) <~ (Tk 2 8)/((3Ha 2 ) where 5 <~ 10~ 5 is the 
density perturbation, k is the wave number and a is the 
scale factor. Using Eqs. (92) and (93), we see that the 
contributions to B a from the vorticity term and the gra- 
dient term are of the same order of magnitude, and we 
obtain 



\B\)^^B^^-^ 



10~ 27 G, (95) 



where we evaluated the amplitude at the horizon scale 
at recombination, H ~ k <~ l/(100Mpc) and a <~ 10~ 3 . 
This is in the range of previous estimates [6, 7, 10, 11], 
as in Eq. (4). 

The anisotropic stress of photons is neglected in the 
present analysis. However, as reported by Ichiki et al. [9, 
11], this is important for magnetogenesis on small scales 
(< 1 Mpc) and in the earlier universe. It is expected that 
the anisotropic stress is important also for the generation 
of vorticity on the same scales and in the same era. We 
will discuss the effect of the anisotropic stresses in future 
work. 



Acknowledgements: 

TK and TS thank Akio Hosoya and Toshio Tcra- 
sawa for their useful comments. RM thanks Chiara 
Caprini, Rob Crittenden, Peter Dunsby, Ruth Durrer, 
Lukas Hollenstein, Christos Tsagas and Ethan Siegel 
for useful discussions. KT is grateful to Kiyotomo 
Ichiki for useful discussion. TK is supported by the 
JSPS under Contract No. 01642. The work of TS was 
supported by Grant-in-Aid for Scientific Research from 
Ministry of Education, Science, Sports and Culture of 
Japan (No. 13135208, No. 14102004, No. 17740136 and 
No. 17340075), the Japan-U.K. and Japan-France Re- 
search Cooperative Program. The work of RM was partly 
supported by PPARC. KT is supported by a Grant-in- 
Aid for JSPS Fellows. 



8 



[1] D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 
(2001) [arXiv:astro-ph/0009061]; 

L. M. Widrow, Rev. Mod. Phys. 74, 775 (2003) 
[arXiv:astro-ph/0207240]; 

M. Giovannini, Int. J. Mod. Phys. D 13, 391 (2004) 
[arXiv: astro-ph/03 12614]; 

M. Giovannini, Class. Quant. Grav. 23, Rl (2006) 
[arXiv:astro-ph/0508544]; 
K. Subramanian, arXiv:astro-ph/0601570; 
R. Durrer, arXiv:astro-ph/0609216. 
[2] See for example, M. S. Turner and L. M. Widrow, Phys. 
Rev. D 37, 2743 (1988); 
B. Ratra, Astrophys. J. 391, LI (1992); 
A. Dolgov, Phys. Rev. D 48, 2499 (1993) [arXiv:hep- 
ph/9301280]; 

M. Gasperini, M. Giovannini and G. Veneziano, Phys. 
Rev. Lett. 75, 3796 (1995) [arXiv:hep-th/9504083]; 
K. Bamba and J. Yokoyama, Phys. Rev. D 69, 043507 
(2004) [arXiv:astro-ph/0310824]; 

K. Enqvist, A. Jokinen and A. Mazumdar, JCAP 0411, 
001 (2004) [arXiv:hep-ph/0404269]; 

K. Bamba and J. Yokoyama, Phys. Rev. D 70, 083508 

(2004) [arXiv:hep-ph/0409237]; 

K. Bamba and M. Sasaki, astro-ph/0611701. 
[3] E. R. Harrison, Mon. Not. R. Astron. Soc. 147, 279 
(1970). 

[4] T. Vachaspati and A. Vilenkin, Phys. Rev. Lett. 67, 1057 
(1991); 

D. N. Vollick, Phys. Rev. D 48, 3585 (1993); 

K. Dimopoulos, Phys. Rev. D 57, 4629 (1998) [arXiv:hep- 

ph/9706513]; 

A. C. Davis and K. Dimopoulos, Phys. Rev. D 72, 043517 

(2005) [arXiv:hep-ph/0505242]. 

[5] A. Rebhan, Astrophys. J. 392, 385 (1992); 

A. Lewis, Phys. Rev. D 70, 043518 (2004) [arXiv:astro- 
ph/0403583]. 

[6] S. Matarrese, S. Mollerach, A. Notari and A. Ri- 
otto, Phys. Rev. D 71, 043502 (2005) [arXiv:astro- 
ph/0410687]. 

[7] R. Gopal and S. Sethi, Mon. Not. Roy. Astron. Soc. 363, 



529 (2005) [arXiv:astro-ph/0411170]. 
[8] K. Takahashi, K. Ichiki, H. Ohno and H. Hanayama, 
Phys. Rev. Lett. 95, 121301 (2005) [arXiv:astro- 
ph/0502283]. 

[9] K. Ichiki, K. Takahashi, H. Ohno, H. Hanayama and 
N. Sugiyama, Science 311, 827 (2006) [arXiv:astro- 
ph/0603631]. 

[10] E. R. Siegel and J. N. Fry, arXiv:astro-ph/0604526. 
[11] K. Ichiki, K. Takahashi, N. Sugiyama, H. Hanayama and 

H. Ohno, arXiv:astro-ph/0701329. 
[12] H. Lesch and M. Chiba, arXiv:astro-ph/9411072; 

C. J. Hogan, arXiv:astro-ph/0005380; 

Z. Berezhiani and A. D. Dolgov, Astropart. Phys. 21, 59 

(2004) [arXiv:astro-ph/0305595]; 

G. Betschart, P. K. S. Dunsby and M. Marklund, Class. 

Quant. Grav. 21, 2115 (2004) [arXiv:gr-qc/0310085]. 
[13] See for example, A. C. Davis, M. Lilley and 

O. Tornkvist, Phys. Rev. D 60, 021301 (1999) 

[arXiv:astro-ph/9904022]. 
[14] G. F. R. Ellis, in Cargese Lectures in Physics, vol. VI, 

ed. E. Schatzman (Gordon & Breach, NY, 1973); 

C. G. Tsagas and J. D. Barrow, Class. Quant. Grav. 14, 

2539 (1997) [arXiv:gr-qc/9704015]; 

C. Caprini, S. Biller and P. G. Ferreira, JCAP 0502, 006 

(2005) [arXiv:hep-ph/0310066]; 

C. G. Tsagas, Class. Quant. Grav. 22, 393 (2005) 

[arXiv:gr-qc/0407080]. 
[15] R. Maartens, T. Gebbie and G. F. R. Ellis, Phys. Rev. 

D 59, 083506 (1999) [arXiv:astro-ph/9808163]. 
[16] P. J. E. Peebles and J. T. Yu, Astrophys. J. 162, 815 

(1970); 

P. J. E. Peebles, The Large Scale Structure of the Uni- 
verse (Princeton University Press, 1980). 
[17] A. D. Dolgov and D. Grasso, Phys. Rev. Lett. 88, 011301 
(2002) [arXiv:astro-ph/0106154]; 

Z. Berezhiani and A. D. Dolgov, Astropart. Phys. 21, 59 
(2004) [arXiv:astro-ph/0305595]. 
[18] K. Takahashi, K. Ichiki and N. Sugiyama, in preparation. 



